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ABSTRACT 

We report Swift Burst Alert Telescope (BAT) observations of the X-ray Flash 
(XRF) XRF 050416A. The fluence ratio between the 15-25 keV and 25-50 keV 
energy bands of this event is 1.5, thus making it the softest gamma-ray burst 
(GRB) observed by BAT so far. The spectrum is well htted by the Band function 
with of 15.oily keV. Assuming the redshift of the host galaxy (z = 0.6535), 
the isotropic- equivalent radiated energy and the peak energy at the GRB 
rest frame of XRF 050416A are not only consistent with the correlation 

found by Amati et ah and extended to XRFs by Sakamoto et ah, but also hll- 
in the gap of this relation around the 30 - 80 keV range of This result 

tightens the validity of the E^^^^ - E^g^ relation from XRFs to GRBs. 

We also hnd that the jet break time estimated using the empirical relation 
between Fp^j^ and the collimation corrected energy E^ is inconsistent with the 
afterglow observation by Swift X-ray Telescope. This could be due to the extra 
external shock emission overlaid around the jet break time or to the non existence 
of a jet break feature for XRF, which might be a further challenging for GRB jet 
emission, models and XRF/GRB unihcation scenarios. 
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1. Introduction 


The observations of X-ray flashes (XRF) are providing important information for under¬ 
standing the nature of Gamma-Ray Bursts (GRB). The detailed studies of XRFs started few 
years ago based on BeppoSAX observations (Heise et al. 2000; Kippen et al. 2002), but X-ray 
rich events had already been detected by the Ginga satellite. Yoshida et al. (1989) reported 
that soft X-ray emission below 10 keV co-exists with y-ray emission of GRBs. About 36% of 
the bright bursts observed by Ginga have F'peak energy, which is the photon energy at which 
the uFiy spectrum peaks, around a few keV and also show large X-ray to y-ray fluence ratios 
(Strohmayer et al. 1998). 

The Wide Field Gameras (WFG) on-board the BeppoSAX satellite observed 17 XRFs 
in hve years (Heise et al. 2000). Kippen et al. (2002) searched for GRBs and XRFs which 
were observed in both WFG and BATSE. The WFG and BATSE joint spectral analysis of 
XRFs shows that their energies are signihcantly lower than those of the BATSE 77°^®)^ 
distribution (Preece et al. 2000). The systematic study of the spectral properties of XRFs 
observed by HETE-2 also supports this result (Sakamoto et al. 2005a). 

The afterglow detection and the redshift measurement from the host galaxy of XRF 
020903, which is one of the softest XRF observed by HETE-2, shows the dramatic progress in 
understanding the nature of XRFs. The prompt emission of XRF 020903 has 77°^!]^ <5.0 keV 
which is two orders of magnitude smaller than that of typical GRBs. The optical transient 
and the host galaxy of XRF 020903 were detected. Further spectroscopic observation of the 
host galaxy suggests that the redshift is 0.25 ± 0.01 (Soderberg et al. 2004). Sakamoto et 
al. (2004) calculated the isotropic-equivalent energy Ei^o and the peak energy at the source 
frame 77®™j^ using the redshift of the host galaxy, and found that XRF 020903 follows an 
extension of the empirical relationship between Ei^o and 77®™j^ found by Amati et al. (2002) 
for GRBs (a.k.a. Amati relation). This result provides the observational evidence that XRFs 
and GRBs form a continuum and are a single phenomenon. 

In this paper, we report the prompt emission properties of XRF 050416A observed by 
Burst Alert Telescope (BAT) on-board the Swift satellite. The X-ray flash, XRF 050416A, 
was detected and localized by the Swift (Gehrels et al. 2004) Burst Alert Telescope (BAT; 
Barthelmy et al. (2005)) at 11:04:44.5 UTG on 2005 April 16 (Sakamoto et al. 2005b,c). 
Swift autonomously slewed to the BAT on-board position, and both Swift X-Ray Tele¬ 
scope (XRT; Burrows et al. (2005)) and UV-Optical Telescope (UVOT; Roming et al. (2005)) 
detected the afterglow (Gusumano et al. (2005) in preparation, Holland et al. (2005) in 
preparation). The afterglow emission of XRF 050416A was also observed by ground obser¬ 
vatories at various wavelengths (Genko et al. 2005a; Anderson et al. 2005; Li et al. 2005; 
Kahharov et al. 2005; Price et al. 2005; Genko et al. 2005b; Soderberg et al. 2005). Genko 
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et al. (2005c) reported that the host galaxy is faint and blue with large amount of the star 
formation and its redshift is z = 0.6535 ± 0.0002. Throughout this paper, the quoted errors 
are the 90% conhdence level and the sky coordinates in J2000 unless we state otherwise in 
the text. 


2. BAT data analysis 

The BAT data analysis was performed using the Swift software package (HEAsoft 
6.0). The background was subtracted using the modulations of the coded aperture (mask¬ 
weighting technique). In this technique, photons with energies higher than 150 keV become 
transparent to the coded mask and these photons are treated as a background. Thus, in this 
mask-weighted technique the effective BAT energy range is from 14 keV to 150 keV. 

Figure 1 shows the energy resolved BAT light curves of XRF 050416A. It is clear that 
the signal of the burst is only visible below 50 keV. The burst signal is composed of two 
peaks. The first peak has a triangular shape with the rise time longer than the decay time. 
When we calculate the spectral lag (Norris et al. 2000) between the 25-50 keV and 15-25 keV 
band, the cross-correlation function lag is —0.066](ogJg second (la error). These temporal 
characteristics are very unusual for the typical GRBs (e.g. Mitrofanov et al. (1996); Norris et 
al. (2000)), thus, it is difficult to understand them in the frame work of the standard internal 
shocks models in which the rise time is always shorter than the decay time and the hard 
emission always exceeds the soft emission (e.g. Piran (1999), Kobayashi et al. (1997)). The 
tgo and tso in the 15-150 keV band are 2.4 and 0.8 seconds, respectively. This fgg belongs to 
the shortest part of the “long GRB” classihcation based on the BATSE duration distribution 
(Paciesas et al. 1999). The fluence ratio between the 15-25 keV band and the 25-50 keV 
band of 1.5 makes this burst one of the softest GRBs observed by BAT so far. The bottom 
panel of hgure 1 shows the count ratio between the 25-50 keV and 15-25 keV bands. The 
spectral softening is clearly visible during the hrst and the second peak. 

As reported by the BAT team^, we applied the energy-dependent systematic error vector 
in the spectral hies before doing any htting procedure. The background subtracted (mask- 
weighted) spectral data were used in the analysis. The XSPEG vll.3.1 software package was 
used for htting the data from 14 keV to 150 keV to the model spectrum. 

Table 1 shows the huences and the peak photon huxes in the various energy bands. These 
huences and peak photon huxes were derived directly from htting the time-averaged and 1-s 


^http://legacy. gsfc.nasa.gov/docs/swift/analysis/bat_digest.html 
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peak spectra respectively assuming the Band function with a = —1. Table 2 summarizes the 
spectral parameters of the BAT time-averaged spectrum^. Figure 2 shows the time-averaged 
spectrum, accumulated over the time interval from —0.5 seconds to 3 seconds since the BAT 
trigger time, was htted with a simple power-law model. The photon index (3 which is much 
steeper than —2 strongly indicates that the BAT observed the higher energy part of the Band 
function (Band et al. 1993). Motivated by this result, and also by the fact that almost all of 
GRB and XRF spectra are well described by the Band function (Preece et al. 2000; Kippen 
et al. 2002), we tried to fit the spectrum with the Band function assuming the low energy 
photon index a to be hxed at —1, which is the typical value for both GRBs (Preece et al. 
2000) and XRFs (Kippen et al. 2002; Sakamoto et al. 2005a). The htting shows a signihcant 
improvement from a simple power-law model to the Band function of of 7.75 for 1 
degree of freedom. To quantify the signihcance of this improvement, we performed 10,000 
spectral simulations assuming our best fit spectral parameters in a simple power-law model, 
and determined in how many cases the Band function £t gives improvements of equal or 
greater than 7.75 over the simple power-law. We found equal or higher improvements in y^ 
in 62 simulated spectra out of 10,000. Thus, the chance probability of having an equal or 
higher Ay^ of 7.75 with the Band function, when the parent distribution is a case of a simple 
power-law model is 0.6%. The observed i7peak energy, i^peak) constrained at I 5 . 6 II 7 

keV, and it confirms the soft nature of this burst. We also applied a constrained Band 
function fit (Sakamoto et al. 2004) to the BAT spectrum to estimate T^peak- The calculated 
T-peak is consistent with the Band function fit of the fixed a to —1: 9.9 keV < F'peak < 20.0 
keV at the 68 % confidence level, 5.1 keV < F'peak < 21.8 keV at the 90% confidence level, 
and i^peak < 23.0 keV at the 99% confidence level. 

The low energy response is crucial for the determination of the spectral parameters of an 
XRF and also, as reported by the BAT team there is a known problem of ~ 15% smaller 
effective area in the Grab spectrum below 20 keV when fitting with a pre-launch response 
matrix. Since the post-launch response matrix which we used in the analysis was applying 
a correction to force the Grab spectrum to fit a canonical model from 14 keV to 150 keV, 
and because we were also also applying the systematic error vectors before performing the 
spectral analysis, the systematic effect of this low energy problem is very limited. However, 
we investigated the spectrum of XRF 050416A ignoring the spectral bins below 20 keV. Even 


^The spectral models which we used throughout this paper are following; a simple power-law model 
(PL): f{E) = KsoiE/SO)^ and the Band function (Band): f{E) = iy3o(if/30)“ exp(— if(2 -f Q;)/ifpeak), 
iiE < {a- d)L;peak/(2 + a) and f{E) = iL3o{(a - d)^^peak/[30(2 + a)]}“-^exp(/3 - a){E/iQf, ii E > 
(a — /3)l^'peak/(2 + a). 


^Section “Corrections to Response” of the BAT Digest (http://legacy.gsfc.nasa.gov/docs/swift/analysis/batjdigest.html) 
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without using spectral bins below 20 keV, the photon index of XRF 050416A is —3.4 ± 0.4, 
much steeper than —2 [a < —2 at the > 99.99% conhdence level). Furthermore, we took 
the ratio of the spectral data of XRF 050416A and the Crab nebula observed at a similar 
incident angle to XRF 050416A. The result is shown in hgure 3. The flattening trend of the 
photon index below 25 keV is also clear in this figure. Thus, we conclude that the deviation 
from a simple power-law model below 25 keV is a real features of the spectrum of XRF 
050416A. 


3. Discussion 

One of the most important discoveries related to XRF 050416A is the conhrmation of 
the F'peak ” -^iso relation (Amati et al. 2002). We calculate the F'peak energy at the GRB 
rest frame, and the isotropic-equivalent energy (1 - 10^ keV at the rest frame), Eiso, 

using the redshift of the host galaxy (z=0.6535). Assuming a = —1, i?pe% and Ei^o of XRF 
050416A are keV and (1.2 ± 0.2) x 10^^ erg, respectively. Figure 4 shows the data 

point of XRF 050416A with the known redshift GRBs of BeppoSAX and HETE — 2 sample 
(Amati 2003; Lamb et al. 2004; Sakamoto et al. 2004). XRF 050416A not only follows the 
-^peak 0^ -^iso relation, but also hlls in the gap of the relation around F'peak of 30 - 80 keV. 
This result tightens the validity of this relation at hve orders of magnitude in E^so and at 
three orders of magnitude in XRF 050416A bridges the gap between XRFs which 

have i?peak of less than 10 keV and GRBs in the F'peak ^ -^iso relation. 

The conhrmation of F'peak ^ -^iso relation from XRFs to GRBs gives us a clear indication 
that XRFs and GRBs form a continuum and are a single phenomenon. There are several jet 
models to explain a unihed picture of XRFs and GRBs. The off-axis jet model (Yamazaki et 
al. 2004; Toma et al. 2005), the structured jet model (Rossi et al. 2002; Zhang & Meszaros 
2002; Zhang et al. 2004), and the variable jet opening angle model (Lamb, Donaghy & 
Graziani 2005) are the most popular models in this aspect. On the other hand, there are 
theoretical models to explain XRFs in the frame work of the internal shock model (Meszaros 
et al. 2002; Mochkovitch et al. 2003) and of the external shock model (Dermer et al. 1999; 
Huang et al. 2002; Dermer and Mitman 2003). The cited jet models and internal/external 
shock models not only explain the existences of XRFs, under certain assumptions, but also, 
in some of their realizations or for some values of their parameters, they can predict the 
F^peak “ -^iso Correlation. 

According to the XRT afterglow observation of XRF 050416A, the decay slope of the 
afterglow emission is ~ —0.9 from 0.015 days to ~ 34.7 days after the GRB trigger without 
any signature of a jet break (Gusumano et al. (2005) in preparation; Nousek et al. (2005)). 
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Using -Upgak and Eiso of XRF 050416A measured by BAT, we can estimate the jet 
break time using the relation between -Upg^k and the jet collimation-corrected energy 
found by Ghirlanda et ah (2004) (Ghirlanda relation). However, there is a debate about 
the assumption of the jet model used by Ghirlanda et ah (2004) to derive the relationship 
between and E^ (Xu 2005; Liang & Zhang 2005). Based on this argument, we use the 
empirical relation between Ei^o, E^eak^ and the jet break time at the rest frame, t|g^, derived 
by Liang & Zhang (2005). Note that there is no assumption of a jet model in the formula 
found by Liang & Zhang (2005), and thus their relation is purely based on observational 
properties. When we use the equation (5) in Liang & Zhang (2005), (i?iso/10^^erg) = 
0.85 X (i?pg^j^/100 keV)^'®^ x (f|g^/l day)“^'^^, the jet break time in the observer’s frame is 
estimated to be ~ 1.5 days after the GRB on-set time. Note that this estimated jet break 
time is consistent with the estimation using the Ghirlanda relation assuming the circum- 
burst density of 3 cm“^. Thus, the estimated jet break time using the empirical E^^^^-Eiso- 
t|gt relation is inconsistent with the null detection of a jet break until more than 34.3 days 
after the trigger by XRT. 

In the off-axis jet model (Yamazaki et al. 2004; Toma et ah 2005), the null detection 
of the jet break in the XRT data of XRF 050416A could be difficult to explain. When we 
assume a bulk Lorentz factor of 100, i^pgak of 300 keV for an on-axis observer, and a jet 
opening angle of 2 degrees, the viewing angle from the jet on-axis is estimated to be ~ 4 
degrees from the observed E^^^^ of 25 keV. According to Granot et al. (2002), when observing 
the jet from an angle two times larger than the jet opening angle, we would expect to see 
a rise in the flux around one day after the burst. It is possible to increase the bulk Lorentz 
factor and to reduce the off-axis viewing angle to achieve the same Doppler factor. However, 
in this case, the afterglow light curve should be close to the on-axis case, thus, we would 
expect to see the jet break around the time we estimated. 

On the other hand, the variable jet opening angle model (Lamb, Donaghy & Graziani 
2005) might work for XRF 050416A if E^ is a constant value. If we assume the values typical 
for GRBs (F^peak = 300 keV and the jet opening angle of 5 degrees), the jet opening angle of 
XRF 050416A is calculated to be 52 degrees because of the inverse relation between i^pgak and 
the jet opening angle in the case in which E^ is a constant. When we used the formulation 
of Sari et al. (1999) applying the estimated jet opening angle, the jet break time will be 64 
days in the case of the circum-burst density of 10 cm“^. Both properties of the low F'peak 
and the null detection of the jet break could be explained in the variable jet opening angle 
model if E^ is constant. However, as Ghirlanda et al. (2004) showed, E^ is not a constant 
parameter, but has a good correlation with When we applied the Ghirlanda relation, 

-^peak 0^ E^’’^, in the variable opening angle model, and re-calculated the jet break time, the 
break time will be 0.7 days assuming the circum-burst density of 10 cm“^. In the variable 
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jet opening angle model, there is no way to explain both the Ghirlanda relation and the null 
detection of the jet break by XRT simultaneously. 

One natural way to explain the non-detection of the jet break feature is that extra 
components are overlaid around a jet break time period. According to the afterglow calcu¬ 
lations in the X-ray band by Zhang et ah (2005), there are several possibilities to hide a jet 
break feature due to some kinds of emission by the external shock. These are the external 
shock emission from 1) the dense clouds surrounding a GRB progenitor (e.g. Lazzati et al. 
(2002)), 2) a moderately relativistic cocoon component of a two-component jet (e.g. Granot 
(2005)) , and 3) a jet with large fluctuations in angular direction (patchy jets; Kumar & 
Piran (2000)). On the other hand, it might be the case that XRFs indeed do not show 
the signature of a jet break in the afterglow. Indeed although the numbers in the sample 
are limited, there is no clear observational indication of a jet break in any XRF afterglow 
light curve so far. If the later case is true, we need to change our view of XRFs completely. 
Thus, the multi-wavelength observations of the XRF afterglows will be crucial to investigate 
whether a jet break feature exists in XRFs or not. 


We would like to thank R. Yamazaki, and D. Q. Lamb for useful comments and dis¬ 
cussions. We would also like to thank the anonymous referee for comments and suggestions 
that materially improved the paper. This research was performed while T.S. held a National 
Research Gouncil Research Associateship Award at NASA Goddard Space Flight Genter. 


REFERENCES 

Amati, L., et al. 2002, A&A, 390, 81 

Amati, L., 2003, GhJAA, Vol. 3, Supplement, pp. 455-460 

Anderson, G., et al. 2005, GGN Girc. 3266 

Band, D.L., et al. 1993, ApJ, 413, 281 

Barthelmy, S.D., et al. 2005, Space Science Review, in press 

Burrows, D., et al. 2005, Space Science Review, in press 

Genko, S.B., et al. 2005, GGN Girc. 3265 

Genko, S.B., et al. 2005, GGN Girc. 3269 

Genko, S.B., et al. 2005, GGN Girc. 3542 



Dermer, C. D., Chiang, J., and Bottcher 1999, ApJ, 513, 656 

Dermer, C. D., and Mitman, K. E. 2003, in ASP Conf. Ser. 312, Third Rome Workshop on 
Gamma-Ray Bursts in the Afterglow Era, ed. M. Feroci et ah (San Francisco: ASP), 
301 

Frail, D.A., et ah 2001, ApJ, 562, L55 
Gehrels, N., et al. 2004, ApJ, 611, 1005 

Granot, J., Panaitescu, A., Kumar, P., Woosley, S.E., ApJJ, 570, 61 
Granot, J., 2005, ApJ, in press (astro-ph/0504254) 

Ghirlanda, G., Ghisellini, G., Lazzati, D., 2004, ApJ, 616, 331 

Heise, J., in’t Zand, J., Kippen, R.M., & Woods, P.M. 2000, in Proc. Second Rome Workshop: 
Gamma-Ray Bursts in the Afterglow Era, Ed. E. Gosta, F. Frontera, & J. Hjorth 
(Berlin: Springer), 16 

Huang, Y. F., Dai, Z. G., and Lu, T. 2002, MNRAS, 332, 735 
Kahharov, B., et al. 2005, GGN Girc. 3274 
Kennea, J.A., et al. 2005, GGN Girc. 3268 

Kippen, R.M., Woods, P.M., Heise, J., in’t Zand, J., Briggs, M.S., & Preece, R.D. 2002, 
in Gamma-Ray Bursts and Afterglow Astronomy, ed. G.R. Ricker & R. Vanderspek 
(Melville: AIP), 244 

Kobayashi, S., Piran, T., Sari, R. 1997, ApJ, 490, 92 
Kumar, P., Piran, T. 2000, ApJ, 532, 286 
Lamb, D.Q. et ah, NewAR, 48, 423 

Lamb, D. Q., Donaghy, T. Q., and Graziani, G. 2005, ApJ, 520, 335 

Lazzati, D., et al. 2002, A&A, 396, L5 

Li, W., et al. 2005, GGN Girc. 3270 

Liang, E., Zhang, B., ApJ in press (astro-ph/0504404) 



- 9 - 

Mochkovitch, R., Daigne, F., Barraud, C., & Atteia, J. L. 2003, in APS Conf. Ser. 312, 
Third Rome Workshop on Gamma-Ray Bursts in the Afterglow Era, ed. M. Feroci et 
ah (San Francisco: ASP), 381 

Meszaros, P., Ramirez-Ruiz, E., Rees, M. J., Zhang, B., ApJ, 578, 812 
Mitrofanov, LG., et ah 1996, ApJ, 459, 570 
Norris, J.P., Marini, G.F., Bonnel, J.T. 2000, ApJ, 534, 248 
Nousek, J.A., et ah 2005 submitted to ApJ (astro-ph/0508332) 

Paciesas, W.S., et ah 1999, ApJS, 122, 465 

Piran, T. 1999, Physics Reports, 314, 575 

Preece, R.D., et al. 2000, ApJS, 126, 19 

Price, P.A., et al. 2005, GGN Girc. 3312 

Roming, P., et ah 2005, Space Science Review 

Rossi, E., Lazzati, D., and Rees, M. J. 2002, MNRAS, 332, 945 

Sakamoto, T., et al. 2004, ApJ, 602, 875 

Sakamoto, T., et al. 2005, ApJ, 629, 311 

Sakamoto, T., et al. 2005, GGN Girc. 3264 

Sakamoto, T., et al. 2005, GGN Girc. 3273 

Sari, R., Piran, T., Halpern, J.P. 1999, ApJ, 519, L17 

Schady, P., et ah 2005, GGN Girc. 3276 

Soderberg, A.M., et ah 2005, GGN Girc. 3318 

Soderberg, A.M., et ah 2004, ApJ, 606, 994 

Strohmayer, T.E., Fenimore, E.E., Murakami, T., Yoshida, A. 1998, ApJ, 500, 873 
Toma, K., Yamazaki, R., Nakamura, T. 2005 submitted to ApJ 
Xu, D. 2005, submitted to ApJ (astro-ph/0504052) 

Yamazaki, R., loka, K., Nakamura, T. 2004, ApJ, 607, L103 



Yoshida, A., et al. 1989, PASJ, 41, 509 
Zhang, B. & Meszaros, P. 2002, ApJ, 571, 876 

Zhang, B., Dai, X., Lloyd-Ronning, N. M., & Meszaros, P. 2004, 601, L119 
Zhang, B., et ah, snbmitted to ApJ (astro-ph/0508321) 


This preprint was prepared with the AAS MRX macros v5.2. 



-11 - 


Table 1: Energy fluences and the peak photon fluxes of XRF 050416A assuming the Band 
function with a = —1 


Energy band 

Energy fluence 

Peak photon flux 

[keV] 

[erg cm“^] 

[ph cm“^ s“^] 

15 - 25 

(1.7 ±0.2) X 10“^ 

2 

^•^-0.3 

25 - 50 

(1.5 ±0.2) X 10“^ 

1.7 ±0.2 

50 - 100 

3.4lJ;° X 10-® 

3.2^°;® X 10-^ 

100 ~ 150 

X 10"® 

X 10"^ 

15 - 150 

(3.5 ±0.3) X 10-^ 

5.0 ±0.5 
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Table 2: The time-averaged spectral parameters of XRF 050416A 


Model 

a f3 

Epeak 

Kso 

xVd-o.f. 



[keV] 

[ph cm“^ s“^ keV“^] 


PL 

-3.1 ±0.2 


(4.3 ±0.3) xlO-2 

50.74 / 57 

PL“ 

-3.4 ±0.4 


(4.7 ±0.5) xlO-2 

43.88 / 53 

Band 

-1 (fixed) < -3.4 

15.6i;f 

3.5tol xlO-^ 

42.99 / 56 


Fitting result using only spectral bins above 20 keV. 



Counts/s/det 



Time since trigger [s] 


Fig. 1.— Light curve of XRF 050416a in five energy bands: 15-25 keV, 25-50 keV, 50-100 
keV, 100-150 keV, and 15-150 keV. The bottom panel shows the hardness ratio between the 
25-50 keV and 15-25 keV band. 
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20 50 

channel energy (keV) 


Fig. 3.— The ratio between the spectral data of XRF 050416A and the Crab nebula. 
The numerator and denumerator of the ratio are the XRF 050416A and the Crab nabula 
spectrum, respectively. The solid line shows the best £t power-law slope of —1.9 derived 
from fitting the data above 25 keV. The bottom panel shows the residuals from this best fit 
power-law slope. The reduced is 7.72 in 20 degree of freedom. 
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Fig. 4.— The isotropic-equivalent energy, Eiso, versus the peak energy at the GRB rest 
frame, for XRF 050416A (red square) and the known redshift GRBs from BeppoSAX 

(circle) and HETE-2 (triangle). The BeppoSAK GRB sample is from Amati et ah (2002) 
and the HETE-2 GRB sample is from Lamb et ah The dotted line is the relation of F^peak 
= 89 (i?iso/10^^erg)°-^ (Amati et ah 2002). 















